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The range of existence and the properties of two essentially different chaotic attractors found in a model of 
nonlinear convection-driven dynamos in rotating spherical shells are investigated. A hysteretic transition 
between these attractors is established as a function of the rotation parameter t. The width of the basins 
of attraction is also estimated. 
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1. Introduction 

Chaotic systems by virtue of their apparently "random" fluc- 
tuations, are expected to frequent, with little bias, all points in 
their phase space over wide parameter ranges. This view im- 
plies that there should not be any abrupt transitions between 
distinguishable chaotic states. Indeed, non-linear transitions 
between attractors are rarely found in turbulent fluids. Con- 
trary to this expectation, a number of examples of discontinuous 
transitions have been observed recently in chaotic and turbu- 
lent fluid systems [1—4]. Bifurcations between distinguishable 
chaotic states appear to be a little less unusual in magnetohy- 
drodynamic flows because of the additional degrees of freedom 
offered by the magnetic field. For instance, various types of 
dipolar, quadrupolar, hemispherical dynamos, and bifurcations 
between them are routinely reported in numerical simulations, 
e.g. [5, 6]. In turn, dipolar dynamos are typically found to 
belong to two distinct regimes - a regime with strong dipolar 
field, and another regime with weaker dipolar component and 
significant multipole contributions e.g. [6, 7]. 

Far more remarkable is a recent finding that two essen- 
tially different chaotic dipolar dynamo solutions may exist at 
identical values of the basic parameters of a generic model 
of convection-driven dynamos in rotating spherical shells [8]. 
Such bistability offers the possibility of a hysteretic transition. 
While hysteresis was established in [8] as a function of three of 
the parameters in the problem - the Rayleigh, the ordinary, and 
the magnetic Prandtl numbers, to be defined below, the depen- 
dence on the last remaining basic parameter, the Coriolis num- 
ber t, was not studied there. This, however, leaves an impor- 
tant gap since the variation in t in this minimal self-consistent 
model of spherical convective dynamos describes the different 



rotation rates characteristic for various cosmic objects. In ad- 
dition, given that current geodynamo simulations are unable to 
achieve geophysically realistic values of t, extrapolation of the 
r dependence is heavily used to compare models and observa- 
tions. As a partial remedy, bistability has been demonstrated for 
two specific values of r in [9]. For these reasons, here we wish 
to investigate the full extent of the coexistence and hysteresis 
in dependence of the rotation parameter r. A number of other 
important questions left open in [8, 9], including the width of 
the basins of attraction of the distinguishable chaotic states and 
the possibility of spontaneous transition will also be discussed 
in the present paper, along with results on essential properties 
of dynamo action such as kinetic, magnetic and cross-helicity 
generation. 

2. Formulation and methods 

2.1. Model 

We employ a minimal model of nonlinear convection-driven 
dynamo process in rotating spherical shells with the aim to re- 
tain general applicability of the results to a variety of cosmic 
objects and to understand fundamental physical mechanisms. 
The model is identical to the one described in [8] but for com- 
pleteness we provide a concise formulation below. We consider 
a spherical fluid shell of thickness d rotating with a constant an- 
gular velocity Q. The existence of a static state is assumed with 
a temperature distribution T$ — To -f3d 2 r 2 /2 and a gravity field 
in the form g - -dyr, where rd is the length of the position vec- 
tor with respect to the center of the sphere. This form of temper- 
ature profile alludes to the possibility that at least a fraction of 
the energy available to planetary dynamos is due to radiogenic 



Accepted for publication in Physica Scripta 



March 2. 2013 



heat release. In addition to d, we use the time d 2 /v, the tem- 
perature v 2 lyad 4 and the magnetic flux density v(jj.g) l ^ 2 /d as 
scales for the dimensionless description of the problem where v 
denotes the kinematic viscosity of the fluid, k its thermal diffu- 
sivity, g its density and /u its magnetic permeability. In common 
with most other simulations of Earth and planetary dynamos 
[10, 11], we assume the Boussinesq approximation implying a 
constant density g except in the gravity term where its tempera- 
ture dependence is taken into account with a = — ( dg/ dT)/g = 
const. The equations of motion for the velocity vector u, the 
heat equation for the deviation from the static temperature 
distribution, and the equation of induction for the magnetic flux 
density B are then given by 

V-h=0, VB = 0, (la) 

(d, + u ■ V)h + rk x u = -Vtt + ®r + V 2 u + B ■ VB, (lb) 

P(d,0 + h ■ V0) = flr • h + V 2 0, (lc) 

V 2 B = P m (d,B + u VB-B V«), (Id) 

where all gradient terms in the equation of motion have been 
combined into Vk. The dimensionless parameters in our for- 
mulation are the Rayleigh number R, the Coriolis number r, the 
Prandtl number P and the magnetic Prandtl number P m , 
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where A is the magnetic diffusivity. Being solenoidal vector 
fields, u and B can be represented uniquely in terms of poloidal 
and toroidal components, 



u = V x (Vv x r) + Vw x r 
B = V x (V/z x r) + Vg x r 



(3a) 
(3b) 



We assume fixed temperatures at r — r t = 2/3 and r — r = 5/3 
and stress-free rather than no-slip boundary conditions in order 
to approach, at least to some extent, the extremely low values 
of viscosity believed to be appropriate to planetary cores [ 12], 



v = dtv = d r (w/r) = = 0. 



(4) 



Two conditions on the poloidal scalar v are required at each 
boundary because the corresponding poloidal equation is ob- 
tained by taking r ■ V x Vx of ( lb) and thus it is of higher order 
as discussed below. For the magnetic field we assume electri- 
cally insulating boundaries at r — r-, and r — r a such that the 
poloidal function h matches the function h ie) which describes 
the potential fields outside the fluid shell, 

g = h- h [e) = d r (h - h {e) ) = at r = r u r a . (5) 

The radius ratio r,/r = 0.4 is slightly larger than that ap- 
propriate for the Earth's liquid core. This is a standard for- 
mulation of the spherical convection-driven dynamo problem 
[10, 11, 13] for which an extensive collection of results already 
exists [5, 6, 14, 15]. The results reported below are not strongly 
model dependent. In particular, dynamos with stress-free and 
with no-slip velocity boundary conditions as well as with dif- 
ferent modes of energy supply are known to have comparable 
energy densities and symmetry properties (see fig. 15 of [ 1 1]). 



2.2. Methods of solution 

Equations of motion for the scalar fields v, w, are obtained 
by taking r ■ V x Vx and r ■ Vx of equation (lb) and equations 
for g and h are obtained by taking r ■ Vx and r ■ of equa- 
tion (Id). These equations are solved numerically by a pseudo- 
spectral method as described in [16] based on expansions of 
all dependent variables in spherical harmonics for the angular 
dependences and in Chebychev polynomials for the radial de- 
pendence. A minimum of 41 collocation points in the radial 
direction and spherical harmonics up to the order 96 have been 
used in all cases reported here which provides adequate resolu- 
tion The dynamo solutions are characterized by their magnetic 
energy densities, 
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M p = - <| V x (V* x r) M t = ^8 X r | z >, 
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-(I V x (V/j x r) | 2 ), M t = -<| Vgxr | 2 >, 



where (■) indicates the average over the fluid shell and h refers 
to the axisymmetric component of h, while h is defined by h = 
h - h. The corresponding kinetic energy densities E p , E t , E p 
and E t are defined analogously with v and w replacing h and g. 
The total magnetic energy density is M = M p + M, + M p + M t , 
and similarly for the total kinetic energy density E. In addition, 
the magnetic energy densities can be divided into those of fields 
that are antisymmetric (axial dipole symmetry) and those that 
are symmetric (axial quadrupole symmetry) with respect to the 
equatorial plane. The former (latter) are described by spherical 
harmonics Y'" with odd (even) Z + m. Other quantities of interest 
are the helicity density of a vector field a, 

H a = a ■ (V X a), 

(known as "kinetic" helicity density when a — u, and "mag- 
netic" helicity density when a — B, respectively [17]), and the 
cross-helicity density 

H x — u • B, 

all of which play important roles in the production of the mag- 
netic field by the chaotic convective flow. 

3. Bistability 

Typical examples of solutions to the dynamo problem out- 
lined above are shown in Figure 1 . The figure presents the main 
magnetic and kinetic energy density components of two dis- 
tinct dynamo cases as functions of time, and illustrates well the 
chaotic nature of the solutions. Apart from the obvious quanti- 
tative difference, an essential qualitative change in the balance 
of magnetic energy components can be observed. The axisym- 
metric poloidal component M p is dominant in the case shown 
in Figure l(a,b) while it has a relatively small contribution in 
the case of Figure l(d,e). This observation is in agreement with 
the claim made in [5-8, 18, 19] that, in general, two regimes 
of dipolar dynamos can be distinguished, namely those with 
M p < M p (denoted by MD, "Mean Dipole" in [8] and below) 
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Figure 1 : (color online) Coexisting distinct chaotic dynamo attractors at iden- 
tical parameter values - a MD (left column (a,b)) and a FD dynamo (right 
column (d,e)) both at R = 1.5 X 10 6 , r = 2 X 10 4 , P = 0.75 and P m = 1.5. 
Panels (a,d) show time series of magnetic dipolar energy densities. The rest of 
the panels show kinetic energy densities in the presence of magnetic field (b,e) 
and after the magnetic field is removed (c,f). The component X p is shown by 
solid black line, while X,, X p , and X, are shown by red, green and blue lines, 
respectively. X stands for either M or E. 
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Figure 2: (color online) (a) Coexistence and hysteresis phenomena shown by 
the time-averaged ratio M p /M p as a function of t. FD and MD dynamos are 
indicated by full red and empty blue circles, respectively. Expected locations 
of the transitions from FD to MD dynamos, and vice versa are indicated by 
dash-dotted lines with arrows pointing down and up, respectively, (b) Com- 
parison of convective properties of FD and MD dynamos measured by the ra- 
tions NuJ D /A r n™ > (violet diamonds), £™/£^ ID (black circles), £f D /£™ (red 
squares), E™ /£™ (green plus signs), E™ /Ef D (blue crosses), all as functions 
of r. Parameter values are P = 0.75, P m = 1.5 and R = (5 — 3 • 
R ■ IQ- S = 7.6, 17, 26, 35, 43, 51, 58, 62 at r = n x 10 4 , n = 1... 
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and those with M p > M p (denoted by FD, "Fluctuating Dipole" 
in [8] and below). The dynamos in Figure l(a,b) and l(d,e) are 
examples of these two types. A convenient measure of the type 
of dynamo is therefore provided by the ratio M p /M p , which we 
use extensively below. 

Far more remarkable is the fact that the two distinct solu- 
tions shown in Figure 1 are obtained at identical parameter val- 
ues and coexist in this case. In fact, this is far from being an 
isolated example. Indeed, varying the value of r we find an ex- 
tended region of coexisting MD and FD dynamos as illustrated 
in Figure 2(a) where the ratio M p /M p is plotted as a function 
of t in the case P = 0.75, P m = 1.5 and R « 4 x R c , where 
R c is the critical Rayleigh number for the onset of convection. 
The transition between the MD and FD dynamos is discontin- 
uous and it is achieved via a hysteresis loop in the following 
sense. When an MD dynamo is used as initial data and the 
Coriolis number r is gradually decreased, solutions remain in 
regime MD until the critical value tfd~ 12500 is reached at 
which point an abrupt jump transition to the FD regime occurs. 
Similarly, when a FD dynamo is used as an initial condition 
and r is gradually increased the reverse transition occurs at the 
critical value tmd~ 39000 as seen in Figure 2(a). The solutions 
plotted in Figure 2 have been typically continued at least up to 4 
magnetic diffusion times. No evidence for a transient nature in 



any case has been found. In fact, in cases outside of the region 
of double attractors it takes typically less than 0.15 magnetic 
diffusion times to switch from the initial conditions used to the 
appropriate unique attractor. 

The results shown in Figure 2(a) are an important comple- 
ment to the findings of [8, 9] where coexistence and hysteresis 
was established as a function of the remaining non-dimensional 
parameters P, P m and R for fixed values of t = 3 x 10 4 , 4 x 10 4 
(see Figures 5 and 10 of [8], and [9], respectively). 

4. Comparison of bistable attractors: Magnetic features 

All solutions included in Figure 2 have a predominantly dipo- 
lar character with the ratio M^' P /M ; „ quantifying dipolarity the 
field visible outside of the spherical shell, in the ranges [0.93, 1] 
for MD dynamos and [0.42,0.70] for FD dynamos, respec- 
tively. Although, the FD dynamos feature an increased contri- 
bution of higher multipoles they are of significant geophysical 
relevance [6, 18, 19]. In this section we wish to discuss some of 
the magnetic properties of MD to FD dynamos in more detail 
than it has been done previously in [8, 9]. 
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Figure 3: (color online) A period of dipolar oscillations of a FD dynamo. The 
plots in the leftmost column are meridional cuts showing lines of constant 
in the left half and of r sin 9dgh in the right half. The plots in the middle column 
show lines of constant B, at r = r„ + 1.3. The plots in the rightmost column are 
meridional cuts showing lines of constant densities of cross-helicity H x , kinetic 
helicity H u and magnetic helicity H B at ip = 0, in the left half, right upper and 
right lower quarters, respectively. The rows correspond to equidistant moments 
separated by At = 0.0252. The parameter values are the same as in Figure 1. 

4.1. Time dependence 

Convection-driven dynamos exhibit chaotic time dependence 
(see e.g. Figure 1) except in simple cases close to the criti- 
cal value of R and for rather large values of P m , such as the 
well-known dynamo benchmark case [20]. There are, however, 
some coherent temporal features that can be distinguished. In 
particular, dynamos in the FD regime are typically oscillatory 
in that nearly periodic changes in amplitude and field structures 
can be observed. As an illustration, Figure 3 shows a period 
of one such dipolar oscillation of the FD dynamo discussed 
in connection to Figure l(d,e). At the beginning of the cycle 
magnetic flux with a "new" polarity is generated near the equa- 
tor of the inner boundary, and subsequently propagates towards 
the poles replacing the flux of "old" polarity, as can be best 
seen in the plots of B v in this figure. The process repeats in a 
quasi-periodic fashion. Kinetic helicity density remains largely 
unaffected while the magnetic and the cross-helicity densities 
participate in the oscillation. On the other hand, dynamos in 




Figure 4: (color online) Time-averaged spatial structures of a MD dynamo. The 
same quantities are plotted as in Figure 3. The parameter values are the same 
as in Figure 1 . 

the MD regime show significantly less variation in time and 
feature spacial structures that remain nearly static and fluctuate 
little with respect to their time averages. This is due to the dom- 
inance of the mean components of the poloidal magnetic field 
characteristic for the MD regime. 

4.2. Spatial structures 

Typical time-averaged spatial structures in the MD regime 
are illustrated in Figure 4 by the example already discussed 
in connection to Figure l(a,b). The dynamo exhibits a nearly 
perfect and relatively large-scale dipolar field best seen in the 
plots of the radial magnetic field and the meridional field lines 
rsinOdgh = const. In particular, two strong zonal magnetic 
flux tubes of By are formed inside the tangent cylinder, near 
the poles, while two tubes of opposite polarity reside on both 
sides of the equator. The kinetic and magnetic helicity densi- 
ties are generated in narrow plumes primarily at the boundary 
of the tangent cylinder, while the cross-helicity density forms 
strong azimuthal tubes of alternating polarity filling the region 
outside of the tangent cylinder. In contrast, the spatial struc- 
tures of FD dynamos, which are exhibited in Figure 3, have 
relatively smaller scales and the evidence of higher multipole 
contributions are clearly visible. At the minimum of the dipolar 
oscillation the zonal magnetic flux tubes near the poles disap- 
pear, and the radial magnetic field shows an excursion towards 
the opposite polarity and an m — 1 structure. The kinetic and 
magnetic helicity densities are again generated near the tangent 
cylinder but appear smaller-scale and more fragmented. 

5. Comparison of bistable attractors: Convective features 

5.1. Non-magnetic convection 

The coexistence and hysteresis phenomena appear to be en- 
tirely magnetic in nature. This is nicely illustrated in Figure 
l(c,f) where the kinetic energy densities of two non-magnetic 
solutions are shown, which are started from initial conditions in 
the corresponding MD and FD dynamos in same figure. Once 
the magnetic field is discarded, the flow quickly equilibrates to 
the same purely convective state. At these values of the pa- 
rameters, the convective state achieved is the well-known state 
of relaxation oscillations [15, 21]. This state of convection is 
remarkable in itself as it has a coherent nearly periodic be- 
haviour in an otherwise chaotic regime. Evidence of its quasi- 
periodicity can be seen both in the time series of Figure l(c,f) as 
well as in the period of relaxation oscillations shown in Figure 
5. 
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Figure 5: (color online) A period of relaxation oscillations seen in the equatorial 
plane of a non-magnetic solution. The top row shows streamlines rd^v = const 
and the bottom row shows color maps of the temperature perturbation & both 
in the equatorial plane. The time step between plots is Af = 0.0448. The same 
case is shown in Figure 1(f). 

Relaxation oscillations are one of a number of states achieved 
by turbulent convection with increase of the Rayleigh num- 
ber R. In this state the differential rotation generated by the 
Reynolds stresses of the convection columns becomes so large 
that it is able to destroy all convective structures by shear- 
ing them off in the azimuthal direction. In the absence of 
convection the differential rotation must decay since there are 
no Reynolds stresses to sustain it. As the shearing action 
of the differential rotation becomes sufficiently weak convec- 
tion columns grow in amplitude again. But as their Reynolds 
stresses regenerate the differential rotation, their amplitude 
quickly peaks and then decays as the shearing action interrupts 
the convection flows. It is surprising how nearly periodically 
this process repeats itself even though every convection episode 
differs from the next one in detail. The period of these relax- 
ation oscillations is primarily governed by the viscous decay of 
the differential rotation. Since our study is based on the viscous 
decay time the value of about 0. 1 was found for the period over 
a wide range of the parameters R, t, and P. 

5.2. Magnetic convection 

The properties of convection in the presence of magnetic 
field differ relatively little between dynamos in the FD and MD 
regime. The time dependence of the respective convective flows 
can be compared in panels (b) and (e) of Figure 1 where kinetic 
energy density components are shown. The time dependence 
is chaotic in both cases, and it may be noted that convection is 
not in the state of relaxation oscillations characteristic for the 
corresponding non-magnetic cases. While the average values 
of the fluctuating components E p and E, are nearly the same 
for the FD and MD case, the main difference is the significant 
increase in differential rotation measured by E,, which in the 
FD dynamo is about 20 times larger than in the MD dynamo. 
Spatial structures of convection for the two dynamos are shown 
in Figure 6. The length scales and the level of irregularity of 
the structures is very similar for the two regimes. The main 
difference appears in the profiles of the zonal flow, which is 




Figure 6: (color online) Spatial structures of convection of a MD and a FD 
dynamos. The leftmost plots show lines of constant tL, in the meridional plane 
on the left, and streamlines rcLv = const, in the equatorial plane on the right. 
The middle plots show lines of constant u, on the spherical surface r = r; + 0.5. 
The rightmost plots show color maps of the temperature perturbation in the 
equatorial plane. The parameter values are the same as in Figure 1. 

nearly constant in z, the direction parallel to the axis of ro- 
tation, in the FD dynamo and strongly dependent on z in the 
MD case. In addition, while the differential rotation near the 
equator is in the prograde direction in the former case, it is in 
the retrograde direction in the latter one. It is well-known that 
the main effect of self-sustained magnetic field on convection 
is to strongly inhibit differential rotation, e.g. [5]. This effect is 
stronger in the case of MD dynamos which are characterized by 
stronger magnetic fields and explains the observed differences. 
Evidence that these differences are typical throughout the MD 
and FD regimes is presented in Figure 2(b) where ratios of var- 
ious kinetic energy components as well as the Nusselt numbers 
are plotted for a number of coexisting attractors as a function of 
t, and it can be seen that their values do not differ significantly 

r • , r -F FD ,-F MD 

from unity except in the case of E t jE t . 

5.3. Mechanism of coexistence 

The observations just made are useful in elaborating the 
mechanism of the coexistence phenomenon. Coexistence is 
the result of two different ways in which the magnetic field 
damps the differential rotation to achieve the transport of the 
same amount of heat. In MD dynamos the differential rotation 
generated by Reynolds stresses of the convection columns is 
eliminated almost entirely by the strong mean magnetic field. 
Because of the strong magnetic field the amplitude of convec- 
tion is also reduced in comparison with the maximum value 
that it reaches in the absence of a magnetic field. In the case of 
FD dynamos the differential rotation is still diminished, but its 
alignment with coaxial cylindrical surfaces is preserved. The 
amplitude of convection is now more strongly fluctuating, but 
is larger on average than in the case of the MD dynamos. In 
this way FD and MD dynamos manage to carry very nearly the 
same heat transport as is evident from Figure 2(b). This heat 
transport by far exceeds the time average of the heat transport 
found in the absence of a magnetic field. 
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Figure 7: (color online) Equilibration of a dynamo started from a small random 
magnetic field seed. Time series of magnetic dipolar energy densities with the 
same color scheme as in Figure 1 are shown. The same parameter values as in 
Figure 1 are used. 
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Figure 8: (color online) Width of the basins of attraction of MD and FD states 
shown by the ratio M p /M p as a function of the continuation parameter a of 
formula (6), in the case R = 1.5 X 10 6 , t = 2 X 10 4 , P = 0.75 and P m = 1.5. 
The cases with a = and a = 1 are the same FD (red full circles) and the MD 
dynamos (blue empty circles) shown in Figure 1 , respectively. 

6. Basins of attraction 

Usually a dynamo equilibrates to a unique state. In contrast, 
for a dynamo inside the region of coexistence, the initial condi- 
tions determine whether a FD or a MD state is approached, as 
discussed in connection with Figure 2(a). Here we attempt to 
estimate the set of initial conditions that leads to equilibration 
to a MD or a FD state. 

6.1. Random initial conditions 

Figure 7 shows a numerical experiment in which random 
magnetic perturbation of small amplitude is applied to fully de- 
veloped convection in the regime of relaxation oscillations, and 
integration in time is continued for the velocity, temperature 
and magnetic fields. After a relatively long transient period of 
magnetic field growth, the dynamo approaches the FD regime. 
In addition, this simulation suggests that the existence of a third 
chaotic attractor along with the FD and the MD states is un- 
likely. 



6.2. Controlled initial conditions 

The vast majority of dynamo solutions published in the liter- 
ature are started from previously equilibrated runs with similar 
parameter values in order to minimize transient times. In the 
region of coexistence a dynamo started from another FD (MD) 
dynamo approaches a FD (MD) state. In order to estimate the 
width of the basins of attraction of the two coexisting states, we 
report in Figure 8 a number of simulations started from initial 
conditions prepared in the form, 

x(r, 6, tp) = ax MD (r, 9,<p) + (l- a)x m (r, 6, <p), (6) 

where a e [0, 1] is a continuation parameter, x represents any 
of the dynamical variables u, B and 0, and the superscripts in- 
dicate equilibrated FD and MD dynamo solutions. When a — 
this is equivalent to initial conditions chosen in the FD regime, 
when a = 1 corresponds to initial conditions chosen in the MD 
regime, and variation of a e (0, 1) allows us to follow a con- 
tinuous path between the two attractors. Figure 8 shows that 
for the test case P = 0.75, R = 1.5 X 10 6 , r = 2 x 10 4 and 
P m = 1.5 the transition between FD and MD regimes occurs at 
a = 0.625 + 0.25. 

We wish to comment that while both Figures 7 and 8 show 
some bias towards the FD regime, this may be due to our in- 
ability to select a test case situated exactly in the middle of the 
coexistence region. 

7. Concluding remarks 

We have considered in this study a minimal self-consistent 
model of dynamo action generated by convection in rotating 
spherical fluid shells. While the relatively thick spherical shells, 
the relatively large values of the Coriolis parameter, and the rel- 
atively low values of the Rayleigh number employed here are 
more appropriate to the problem of geo- and planetary mag- 
netism, the model is generic and may be used to understand 
solar and stellar magnetism. For instance, the periodic rever- 
sals characteristic for FD dynamos are reminiscent of the 11- 
year Solar cycle. In particular, we have been concerned here 
with the possibility of coexistence of two nonlinear attractors in 
the fully-developed chaotic dynamo regime, and the hysteretic 
transition between them. These phenomena have been noted 
previously in [8, 9] and discussed there in some detail. The 
following is a summary of the important points made here. 

(a) We have established the coexistence of two nonlinear at- 

tractors, denoted by MD and FD dynamos in the above, 
over a significantly large interval of values of the Coriolis 
number t e (12500, 39000). The transition between them 
takes the form of a hysteresis loop. These results fill a gap 
left in [8, 9] and demonstrate that coexistence occurs as a 
function of all basic parameters in the model. 

(b) We have discussed in detail the contrasting properties char- 
acterizing dynamos in the MD and FD regimes, including 
differences in temporal behaviour and spatial structures of 
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both magnetic field and convection. We include new re- 
sults on quantities important in mean-field dynamo theo- 
ries of magnetic field generation such as the kinetic, mag- 
netic and cross-helicity density profiles in time-averages 
as well during oscillations and reversals. 

(c) We have investigated the question of the width of the basins 
of attraction of the coexisting chaotic states. 

The coexistence of two distinct turbulent attractors is also a 
phenomenon of general interest as it is relatively rare in fluid 
dynamics and magnetohydrodynamics. Finally, the range of 
values of the Coriolis number t where we have found coexis- 
tence constitutes a rather large subinterval of the range currently 
accessible by numerical simulations. This requires extra care 
when numerical results are interpreted. 
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